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Abstract 

We propose an exact method for locating the zeros of the Jost function for 
analytic potentials in the complex momentum-plane. We further extend the 
method to the complex angular-momentum plane to provide the Regge tra- 
jectories. It is shown, by using several examples, that highly accurate results 
for extremely wide as well as for extremely narrow resonances with or without 
the presence of the Coulomb interaction can be obtained. 

PACS numbers: 03.65.Nk, 03.65. Ge, 21.45. +v 

I. INTRODUCTION 

Much effort has been devoted in the past to develop methods to calculate the energies 
and widths of resonances in the potential scattering theory. A comprehensive survey of this 
subject can be found in Ref. Jl|. These methods can be divided into two categories. The 
first one is based on methods traditionally employed for real energies where one can locate 
the position of relatively narrow resonances with a sufficiently high accuracy, but has many 
difficulties in determining their widths and usually fails for broad resonances. In the second 
one the calculations are performed at complex energies, and therefore the widths and reso- 
nant energies are obtained simultaneously. 



The complex methods have the advantage that the calculations are based on a rigorous 
definition of resonances, namely, as singularities of the S-matrix. Thus, in addition to the 
information about the resonances themselves, they provide us with information about the 
analytic properties of the S-matrix and the off-shell properties of the underlying interaction. 
However, the existing complex-energy methods are much more complicated than those of 
the first group and require sophisticated mathematical and computational methods to han- 
dle them. 



Usually, the complex methods are referred to as the "direct calculation approaches" , but 
very often with the quotation marks J2| because most of them are based on an expansion 
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of the resonant wave function in terms of square integrable functions and the subsequent 
determination of the expansion coefficients either by a diagonalization of the non-Hermitian 
Hamiltonian or by a variational procedure. 

The method we present here formally belongs to the second group i.e., to the complex- 
energy methods. It is based on exact differential equations for functions closely related to 
the Jost solutions and which coincide with the Jost functions at large distances Un- 
like the existing methods, it is simple in applications and although it exploits the idea of 
the complex rotation of the coordinates, it is different from the traditionally used complex 
dilation methods in that it does not employ expansion or variational procedures. Instead, 
the Jost function at a complex energy is obtained directly from exact equations equivalent 
to the Schrodinger equation. 

To demonstrate the effectiveness and accuracy of our approach, we performed calcula- 
tions using potentials previously employed by other authors We not only located the 
resonances cited by them but also find sequences of Jost function zeros corresponding to 
broad and extremely narrow resonances which were not considered or missed in the afore- 
mentioned references. 

In addition to the location of zeros in the complex fc-plane, the present method en- 
ables us to locate the zeros of Jost functions in the complex angular momentum plane. We 
demonstrated this by locating the zeros, known as Regge poles, representing resonances in 
the complex £-plane. 

The paper is organized as follows. In Sec. II our formalism is presented, while Sec. Ill is 
devoted to the boundary conditions. In Sec. IV the method is applied to several examples 
and the results obtained are discussed and compared with those obtained by other methods. 



II. BASIC EQUATIONS 



Consider the one-channel problem of two, generally charged, particles. Appart from the 
Coulomb force, we assume that these particles interact via a central potential V(r) having 
the properties 



lim r 2 V(r) = 



and 



lim rV(r) = . 

r^oo 

The radial Schrodinger equation reads {% = 1) 

d r + 2mk 2 - £(£ + 1) /r 2 - 2r)k/r] $ e (k, r) = V(r) ® e (k, r) 



(1) 



(2) 



(3) 



The regular solutions <£v(fc,r), for any complex k ^ and complex I in the half-plane 
Re£ > —1/2, are defined by the boundary condition 
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lim [®i(k,r)/F e (r),kr)} 



(4) 



where Fe(rj, kr) is the regular Coulomb function (1C . 



In contrast to $e(k,r), the physical solutions (bound, scattering, and Siegert states) are 
defined by their behaviour at large distances. However, they all are regular at r = and 
thus proportional to $>i(k, r). Therefore, once the function $^(A;,r) is determined at all 
complex momenta k, it represents, in a most general form, all solutions of physical interest 
of Eq. (H|) since any physical solution at specific values of k can be obtained from it merely 
by multiplying r) by the proper normalization constant. 



At large distances $>e(k,r) can be written as a linear combination 



$t(k,r) 



Hi + \ V , kr)f;( V , k*) + H ( e -\ V , kr)f e ( V , k) 



Hi 



(5) 



The functions H i (r],z) are defined in terms of the regular Fi(r),z) and irregular Ge(rj } z) 

(6) 



Coulomb functions flQjlllJlql ), 

H { e ±) ( V ,z)=F e (r ] ,z)T^G e ( V ,z 
and have the asymptotic behaviour, 



^ (±) (?7, z) — ► =Fzexp{±z [z - hi 2,2 - £ir/2 + argr(£+ 1 + irj)]} 



\z\—*oo 



(7) 



For neutral particles or high energies where the Sommerfeld parameter 77 — > 0, the Coulomb 



functions reduce to the Riccati-Bessel, Riccati-Neumann, and Riccati-Hankel functions ||10 
i.e., 



The momentum-dependent coefficients in the above linear combination, Eq. (|5D, are the 
Jost functions which are analytical for all complex k of physical interest and for Re£ > —1/2. 
The last restriction on i stems from the fact that at Re£ = —1/2 the role of Fe(rj,z) and 
G e(r], z), of being regular and irregular at the origin, is interchanged jl I 



For integer (physical) £ the Jost function has zeros at a discrete sequence of points k 0i , 
i = 1, 2, 3, • • • , situated either on the imaginary axis of the fc-plane or under the real axis. 
At these points only the first term in the asymptotic form (|) survives, corresponding to 
either a bound (Im k > 0) or a Siegert (Im k < 0) state behaviour for large r. 

On the other hand, for real values of k 2 (physical energies), the function fe(i], k) can have 
zeros at complex i which move, with increasing k 2 , along the so-called Regge trajectories 
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which define families of bound and resonant states [fLl]] . Therefore, when the Jost function is 
known at all complex values of k and at all permissible values of £, it contains all information 
and characteristics of the spectrum of the underlying physical system. 



In Ref. we proposed a method for a direct calculation of the Jost function for any com- 
plex momentum of physical interest. In this approach we use a combination of the variable- 
constant [TjJ and the complex coordinate-rotation (TT3| methods to solve the Schrodinger 
equation (^|) in an efficient and accurate way without resorting to any approximation, ex- 
pansion, or to variational (stabilization) procedures. For this, we perform a complex rotation 
of the coordinate r 



x exp(zt 



x > 0, 9 G [0, 9 r , 



9„ 



< tt/2 



in the Schrodinger equation (Q) and look for a solution in the form 



$i(k,r) 



1 r 



Hp\r],kr) J^faJfe.z.G) + H\ >(r],kr) Ff'(r),k,x,9) 



:(+) 



(8) 



(9) 



where J-g{rj,k,x,9) are new unknown functions (variable constants) which at large x be- 
come, according to Eq. (|), true constants. In this way the initial value problem, defined 
by Eqs. (0) and (HI), reduces to the following first-order coupled differential equations 



d x ^ + \ V ,k,x,9) 



2ik 



H^\r],kr) V(r) 



H { e +) (77, kr) T^ 1 (77, k, x, 9) + Hf (77, kr) Tf> (r), k, x 



r(+) 



(-)/ 



:(-)/ 



(10) 



d x ^f\ V ,k,x,9) = -—Hf\ V ,kr) V(r) 



?: } %k,x,9) g 

x "ifj"%,Ax) ^ + \ri,k,x,9) + H^\ri,kr) jj~\ri,k,x, 
with the simple boundary conditions 

which follow immediately from Eqs. (13), (0), and (M). 



(12) 



In Ref. it was proved that if the potential obeys the conditions (HJ) and (0), with 
complex r defined by @, then for all complex k situated above the dividing line \k \ exp(—i9) 
and for x — > 00, the function T\ >(tj, k, x, 9) has a ^-independent limit which coincides with 
the Jost function, i.e., 



lim T\ \t), k, x, 9) = f e (r], k) 



(13) 



while the function J 7 ^ (r), k, x, 9) converges to fl(T],k*) at all spectral points k 0ij i 
1, 2, 3, • ■ -, for which f t (r), k 0i ) = 0. 
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The proof was based on the asymptotic behaviour of the functions H\ (r/, kr) at large 
r, and can be generalized to include the complex angular momentum i as well. This gen- 
eralization is straightforward since the angular momentum appears only in the phase factor 
of the asymptotic form, Eq. ([?]), and therefore for a complex t the functions T\ ; (77, k, x, 9) 
at large x have the same behaviour and the asymptotic relation ( |IBD is also valid. 

Thus, the procedure of calculating the Jost function is very simple since for any fixed 
pair of k and I we only need to solve the system of first-order differential equations (|10D 
and (PH) from x = to some x max where t\ \rj,k,x,9) attains a constant value (usually, 
x m ax is the range of the potential V). According to Eq. flT3|), this constant is just the Jost 
function fe(f], k) we are looking for. Simultaneously, as a bonus, we obtain the exact wave 
function in the form of Eq. (^) having the correct asymptotic behaviour. Depending on the 
choice of the momentum k, it can be a bound, scattering, or a Siegert state wave function 
(rotated when 9 > 0). 

The resonances in a specific region of complex k can be easily located by taking the rota- 
tion angle 9 large enough to cover this region and then search for zeros of the Jost function. 
The zeros in the complex £-plane can be similarly located with 9 = 0. 

We emphasize that this method is exact. Although we employ the complex rotation, we 
do not need any stabilization procedure. This has been demonstrated in Ref. M where we 
employed an analytically solvable model and showed that Eqs. ([[(]) and (pTTJ ) give at least 
5— digit accuracy for the Jost function in a wide area of complex k despite the fact that the 
simplest Runge-Kutta method of integration was used. 



III. BOUNDARY CONDITIONS 
A. Short distances 



Formally, we have to start the integration of Eqs. (|T0|) and (JTT]) from x = 0. However, 
for i 7^ the functions (j], kr) are irregular, i.e., at the origin they behave as [14 



where 



TT^fr, hr\ — > T% I — I + 1° (Vkrlnkr) , for i = 

n e {V , Kr) ^ + ^ Ci ^ y 2 j ^\ Q^ikr) 1 -"), for 



In Ref. it was shown that the corresponding singularities at x = in the above differential 
equations (|T0D and ( pTTj ) are integrable when the condition ([!]) is fulfilled. Thus, there is no 
problem from a formal point of view. However, in practical calculations we cannot start 
from x = and therefore we have to shift the initial point to some small value x m i n . Thus, 
to implement the boundary conditions, we need to know ^^(rj, k, x min , 9). 
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There are several ways to circumvent this problem. One of them consists in transforming 
the differential equations, Eqs. (0) and (|TT1), into an equivalent pair of integral Volterra- 
type equations, viz., 

^ {±) (r/, k,x,6) = 1 ± ^ £ HF\ri,kx , e*)V(xre*)$ l (k,x'e*)<y , (14) 

where $>g(k, r) is defined by Eq. (|9|). We can solve these integral equations in the interval 
[0, x m i n ] iteratively as follows: 



J~i ^ \f], k, x m i n , 6) — 1 , (15) 
^ ±)(1) fa,Mm*i,0) = 1±C [ Xmm HF\7 i ,kxe i0 )V(xe i0 )F i (ri,kxe ie )dx , 



rt ±HN \v, k, x mm , 0) = l±^ H^\ V , kxe ie )V(xe ie 



and then integrate the differential equations starting from the value of (rj, k, x m i n , 9). 

For small values of x m i n the above iteration procedure converges very fast. Moreover, in 
implementing the method we found that if Re£ is small (~ 1), then a surprisingly high ac- 
curacy (better than 7— digits) can be achieved even with the lowest order iteration, Eq. (p~5|) . 
For higher values of Re £, however, the use of these simple boundary conditions could result 
in numerical instabilities. This is due to the ansatz (||) which is suitable for large distances, 
but is not good in the vicinity of r = 0. Indeed, near this point the function $>g(k, r), by its 
definition, is regular and therefore the singularities of H^ + \r],kr) and H^~\r],kr) are can- 
celled. This is secured by the boundary conditions (|12|)- In numerical calculations, however, 
the cancellation of singularities is always a precarious procedure and a source of possible 
numerical errors. These errors increase with increasing Re£ since in this case H\ (rj, kr) is 
more singular. Therefore, the greater Re I is the further the point x m i n must be shifted from 
the origin in order to avoid cancellation errors. This shift in turn, requires more iterations 
of Eq. fli~4]) to obtain the boundary values J-£ (tj, k, x min , 6) to a required accuracy. 

Another way to handle the boundary condition problem is to replace at short distances 
the ansatz (^) by a more suitable one. Indeed, this ansatz was motivated by the variable 
constant method ||12|| , i.e. we looked for a solution of Eq. (^) in the form of a linear com- 
bination of its two independent solutions Hg and Hjf~^ corresponding to V(r) = 0. When 
the potential is taken into account, the coefficients of this combination are r-dependent and 
obey the Eqs. (0) and ([□]). Thus, instead of H^^^kr) we can choose another pair of 
linearly independent solutions, namely, F^r), kr) and Ge(t], kr), and the new ansatz reads 

$ e (k, r) = F e (r], kr)A e (ri, k, x, 9) + Gfa kr)B e (r], k, x, 9) . (16) 

Since (|9|) and ([16]) are merely different representations of the same function, we have 
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T) ± \r l ,k ) x,9) = A l {r l) k,x,9)±iB l {7 1 ,k,x ) 9) ) (17) 
and the equations for the functions Ag(r), k, x, 9) and Bg(rj, k, x, 9) are 
d x A e (r],k,x,9) = 

ie 

— G e (ri, kr) V{r) [F e ( V , kr) A e (r,, k, x, 9) + G e (r,, kr) B £ {r], k, x, 9)] , (18) 

d x B e (r],k,x,9) = 

ie 

- ~i-Fl{ri, kr) V{r) [F e ( V , kr) A t {r], k, x, 9) + G t {r), kr) B e (rj, k, x, 9)} . (19) 
The corresponding boundary conditions, 

Ag(r), k,0,0) = 1, B t (n, k, 0, 9) = , (20) 
follow immediately from (|T^) and (JT7|). 

In other words, we have two equivalent pairs of equations, Eqs. (]T0f) and f[TTl) and (|18|) 
and (0), defining the same function <&e{k, r) in its two different representations @ and (|16|). 
Computationally it is easier to start the integration of equations ( |18|) and (|l9f) at x m j n and 
continue it up to some intermediate point Xi n t (not necessary small), and then to integrate 
the equations ([10]) and ([TT|) from a;^ to x ma;c where J-j (r), k, x max , 9) coincides with the 
Jost function. 



Similarly to equations (|10|) and ($A\) the differential equations for Ai(r],k,x,9) and 
Bi(rj, k, x, 9) can be transformed into integral Volterra-type equations, 

A e (r], k,x,9) = l + ^- [ X GArj, kx'e w ) V(x'e i9 )$ e (k, x'e i9 )dx' , (21) 
k Jo 

BArj, k,x,9)= - e — f X FArj, kx'e ie )V{x'e ie )<5> e {k, x'e w )dx' , (22) 
k Jo 

where &e(k, r) is defined by Eq. flip]). Iterations of these integral equations can also be used 
to obtain corrections, if necessary, to the simplest form of the boundary conditions, namely, 

A e (ri,k,x min ,9) = I , B e (r), k,x min ,9) = 0. (23) 



B. Large distances 

The behaviour of <lv(A;,r) at large distance, is determined by the functions H^^^kr). 
Therefore, the correct asymptotic form is automatically secured. 

Indeed, suppose we have found on the positive imaginary axis of the fc-plane a value ko 
for which Fg(j], ko, x max , 0) = (when Im k > we can always put 9 = 0), i.e., we located 
a zero of the Jost function corresponding to a bound state. The physical bound state wave 
function is then given by, 

if b r nd (k ,r)=Af^(k ,r), 
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and differs from &e{ko, r) only by a normalization factor TV which can, in principle, be found 



along with $^(&o,r) in terms of the Jost function and its derivative [15,16], or simply from 



the normalization integral. At large r only the first term of Eq. (^]) survives, i.e., 

$/(fco,r) ,— > ^ + \v,k ,x maxi 0)H { e + \ V ,k r). (24) 

r|— >oo z 

Obviously, in this expression the exponentially decaying tail of the bound state wave func- 
tion is presented in an exact form. 

For scattering states (real positive k), the asymptotic form of $e(k, r) contains both 
terms of Eq. (|9]), 

$ e (k,r) — > I \Hi + \rj, kr) jf~\r}, k, x max , 0) + H^\v, kr) ^ H (?7, k, x max , 0)1 , (25) 

\r\— >oo ^ L J 

where the functions (j], kr) represent the incoming and outgoing spherical waves (again 
in the exact form). The scattering wave function, y?£ C jf*(?"), differs from &i(k,r) only by a 
constant factor, viz., 



i {Vi kj x max , 0) 



HF'fa kr) (.V, k, x, 0) + Hf } (rj, kr) (rj, k, x, 0)J , 

(26) 

where we assumed that the scattering states \ty s S att ) are normalized according to 

^scatt^spattj = - k) 

and expanded in partial waves as follows: 

' im 

For the Siegert states, corresponding to zeros of T\ (r),k,x max ,9) in the lower half of 
the fc-plane, we have the same kind of asymptotic behaviour as in ^(fco, r) given by Eq. 



(24]) but in this case 9 > 0, and therefore such states can be treated and normalized similarly 



to bound states 0. 

In summary, the representation (0) secures the proper behaviour of $>g(k,r) at short 
distances, while the representation (^]) guarantees its correct behaviour at large r. The use 
of these representations enables us to achieve high accuracy in the solution of Eq. @ at all 
complex values of k. 



IV. NUMERICAL EXAMPLES 

In order to demonstrate the ability and accuracy of the proposed method, we chose two 
simple potentials previously used in Refs. |6|-|9|]. This choice is further motivated by the 
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richness of the spectra generated by these potentials and by their simple form. And, as we 
found, their spectra include very wide as well as extremely narrow resonances which are 
difficult to locate with most of the existing methods. Thus they are ideally suited as testing 
cases. 



In atomic units [17], these potentials have the following form 



and 



Vi(r) = 7.5r 2 exp(— r) H — 



V 2 (r) = 5exp [-0.25(r - 3.5) 2 ] - 8exp (-0.2r 



The reader not accustomed to the atomic units, may assume that the above potentials are 
given in MeV and the distances in fm. In such a case h 2 /2m = 1/2 MeVfm 2 while the 
Sommerfeld parameter is given by rj = z/k. Then the numerical values of the resonance 
energies are the same and independent of the unit used (MeV or atomic units a.u.). In what 
follows, in order to avoid possible misunderstanding, we will use the MeV-fm units. We note 
that V\(r) is a good case to test the ability of the method to deal with interactions having 
a Coulomb tail. Similarly to Ref. || we assumed that the Coulomb part is attractive with 
z — — 1. In order to compare our results with those given in Refs. ||-§§, we also performed 
calculations with z = 0. 



To locate zeros of the Jost function in the complex fc-plane as well as in the complex 
£-plane, we searched for minima of its modulus, \^~\r], k, x maxi 9)\, considered as a func- 
tion of two variables, either Refc and Imk or Re^ and lm£. This is based on the so called 



maximum modulus principle for a complex-valued function [18| . According to this principle, 
when a function f(z) is holomorphic and not constant in a region D of the complex plane, 
\f{z)\ can never attain its maximum in the interior of D but only on the boundary of D. 
Therefore the minima of \f{z)\ must coincide with the zeros of f{z). Indeed, assuming that 
\f(z)\ has a minimum at z = zq inside the area D and f{zo) ^ then around the point z 
the function l/f(z) is holomorphic and has a maximum inside an area, which contradicts 
the maximum modulus principle. Thus, if a minimum of \f{z) \ is found, it must be the zero 
af/(z). 



The located zeros of the Jost function in the fourth quadrant of the complex fc-plane 
for the potential V\(r) for the i = partial wave, and the corresponding resonance energies 
and widths are presented in Table I together with the results of Refs. 0-0. Our calcula- 
tions have been performed with the simplest boundary condition (^). This was sufficient 
to achieve an accuracy of at least 9 digits. This is checked by changing the rotation angle 9 
since the Jost function must be ^-independent. We note that only few S'-wave resonances 
are presented although many more were located. The reason is that this potential has been 
used earlier by several authors (but only for the case t = 0) and thus we employed it in 
order to compare our results with those of Refs. 0-0. The sequence of the S'-wave resonant 
zeros continues downwards in the complex fc-plane. This behaviour can be seen in Fig. 1 
where these zeros are plotted. A similar behaviour was also found for resonances of higher 
partial waves. In these cases, however, the first few zeros of the sequences are closer to the 
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real axis (which means that they have a smaller width). The same trends were found for 
the resonances of the potential V^(r) which generates a richer spectrum. 

Comparison with the other calculations for the potential V\(r) shows (see Table I) that 
only the "direct" (dilatation) method of Ref. || gives an accuracy which is comparable with 
ours. In that reference, the 5 digit accuracy was achieved by using over 40 exponentially 
decreasing functions in the expansion of the rotated Siegert state. In contrast a 9 digits 
accuracy is achieved by our method without any exertion and, if neccessary, can easily be 
improved. Such an improvement is of crucial importance when one deals with extremely 
narrow resonances. When the width of a resonance is 7 orders of magnitude less than its 
energy, one needs at least a 7— digit accuracy to be able to discern it. 

Another extreme situation is the case of very broad resonances, i.e., when the Jost func- 
tion zeros are situated far below the real axis. As can be seen in Table I, even for the first 
resonance (E = 3.426390331 MeV, T = 0.025548962 MeV) only a third digit accuracy in 
the width has been achieved by the "real-energy methods" of Refs. J7||| and by the "semi- 
complex method" of Ref. || while the next, moderately broad, resonance (Eq = 4.834806841 
MeV, T = 2.235753338 MeV), is already beyond their resolution. 

The "real-energy methods" of Refs. |7||| consider eigenenergies of a system enclosed in 
a box. These eigenenergies are moving together with the change of the radius of the box, 
generating the so-called quasi-crossings at resonance energies. The width of a resonance is 
determined by the breadth of the quasi-crossing. Such an approach exploits the physical 
idea that the resonant states, are only slightly affected by variations of the radius of the box, 
in contrast to the spurious states that emerge in the box. It seems, however, that this idea 
is not suitable for broad resonant states. Indeed, in Fig. 1 of Ref. and Fig. 1 of Ref. |§ , 
presenting the box eigenenergies, no traces of the the second quasi-crossing of the 5-wave 
resonance at E = 4.834806841 MeV or of the third at E = 5.277279780 MeV can be found. 



It is claimed that broad resonances are unimportant and thus the inability of a method 
to describe them is a minor drawback. However, in certain physical systems, broad reso- 
nances play a significant role. An example is the Sii(1535) resonance of the interaction of 
the ?7-meson with a nucleon which lies 48 MeV above the threshold while its width is 150 
MeV [|19|| . Nevertheless it prescribes the dynamics of the 7/-meson interaction with nuclei. In 



particular, due to this resonanance, certain 77-nucleus systems can have quasi-bound states 



One of the advantages of the exact method presented in this work, is that bound, scat- 
tering, and resonant states can be treated in a uniform way regardless of their width. All 
spectral points can be located with the same accuracy irrespective of their location on the 
complex /c-plane. This is exemplified by the spectral analysis of the potential V^r). Se- 
quences of the bound and resonant states generated by this potential are given in Tables II 
and III. The S- and P-wave sequences are also shown in Fig. 2. 



The S"-wave states, generated by the potential V 2 (r), were previously considered in Ref. 



10 



0, where only the first two resonances were found. As can be seen in Table II, the spherical- 
box method of Ref. j7[ provides the position of the first (narrow) resonance fairly well. How- 
ever, the width is obtained only to two correct digits. For the second resonance, which is 
broader than the first one, the accuracy for the energy evaluation is down to 3 digits and for 
the width just to 1 digit. The other resonances of the .S-wave sequence, presented in Table 
II, were not obtained with the box-method. 

The accuracy of the present method is determined by the accuracy of solution of the 
differential equations. By choosing the tollerance to be small enough we can locate prac- 
tically all resonances. Two remarkable examples of extremely narrow resonances, gen- 
erated by the potential V^r), are the one found for the P-wave at Eq = 0.807634844 
MeV (r = 0.110 x 10~ 6 MeV) and the other for the F-wave at E = 1.009031953 MeV 
(r = 0.46 x 10" 7 MeV) and are parts of the sequences presented in Table III. If the attrac- 
tion of the potential is slightly increased, these extremely narrow resonances become bound 
states in the P and F partial waves. 



Some of the Jost-function zeros found in the complex £-plane are given in Tables IV 
and V. The corresponding Regge trajectories are depicted in Fig. 3. Each trajectory begins 
from an S-wave spectral state. The first one begins from the ground state and passes via 
the lowest states of each partial wave. When the energy is negative, the trajectory lies on 
the real axis, while at positive energies it gradually goes upwards. We note that the width 
of a narrow resonance can, in principle, be found via lm£ which corresponds to an integer 
value of Re £. However, the relation between T and Im I involves also the derivative of the 
trajectory [[16 . 



The Regge trajectories, can provide us with a useful information on the spectrum of 
the physical system. They combine bound and resonant states in families and therefore in 
calculating such trajectories, we can find out at which energies and in which partial waves 
resonances must exist. Furthermore, the Regge trajectories give us a general insight into 
the spectrum of the Hamiltonian under consideration. For example, we can state that the 
F-wave resonance at E = 1.009031953 MeV is the narrowest one for the potential V^r). 
This follows from the fact that this resonance belongs to the lowest Regge trajectory and is 
the first on it. 



In conclusion, the method we present in this article, enables us to locate the zeros of the 
Jost function in the complex fc-plane and to calculate the Regge trajectories in the complex 
£-plane in a simple, efficient, and extremely accurate way. To the best of our knowledge, no 
other method achieved such a performance in the past. Since it is based on exact equations, 
the method simultaneously provide us with the corresponding wave functions of the spectral 
states. 



As a final note we mention that the method can be easily generalized to treat the coupled 
channel problems having the same angular momentum. In such a case we replace the 
potential and the functions ^ (?7, k, x, 9) by matrices and the spectrum is then defined 
by the zeros of the Jost-matrix determinant. Channels of different angular momenta can 
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also be treated in the same manner, but this requires a more elaborate treatment of the 
boundary conditions at x = since some off-diagonal elements of the matrices J^ (±) for 
different £ may be singular at the origin. Nonlocal potentials can also be considered within 
the proposed approach. Work on all these generalizations is under way. 
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TABLES 



TABLE I. The zeros ko of the Jost function in the complex fc-plane and the energies Eq and 



widths L of the S-wave resonances for the potential V\(r) with [z = —1) and without (z = 0) the 
Coulomb term. 



z 


k 


(fm- 1 ) 


E Q (MeV) 


L (MeV) 


Ref. 





2.617786172 


- i0.004879876 


3.426390331 


0.025548962 


this work 








3.42639 


0.025549 


1 








3.4257 


0.0256 


i 








3.426 


0.0256 


I 








3.426 


0.0258 


i 





3.130042444 


- i0.357144253 


4.834806841 


2.235753338 


this work 





3.398392393 


- i0.997251873 


5.277279780 


6.778106356 


this work 


-1 


1.887074210 


- i0.000025362 


1.780524536 


0.000095719 


this work 








1.7805 


0.0000958 




-1 


2.871167766 


- i0.201530270 


4.101494946 


1.157254428 


this work 


-1 


3.169186525 


- i0.846652839 


4.663461099 


5.366401539 


this work 



TABLE II. The zeros of the Jost function in the complex £;-plane and the energies Eq and 
widths L of the 5-wave bound and resonant states for the potential V^r). 



i 


ko (fm- 1 ) 


Eq (MeV) 


L (MeV) 


Ref. 





i3.023634507 


-4.571182814 





this work 





il.329872758 


-0.884280776 





this work 





2.122442334 - i0.000027859 


2.252380731 


0.000118256 


this work 






2.25237 


0.0001196 







3.000600515 - iO.041316851 


4.500948186 


0.247950731 


this work 






4.50 


0.28 


% 





3.485234669 - i0.360967988 


6.008281406 


2.516116297 


this work 





3.974580284 - i0.788297975 


7.587937367 


6.266307179 


this work 





4.454733926 - i\. 227097054 


9.169443586 


10.932781876 


this work 





4.922800349 - i\. 664647328 


10.731456273 


16.389452891 


this work 





5.378677040 - ^2.097566794 


12.265190122 


22.564268707 


this work 
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TABLE III. The zeros ko of the Jost function in the complex /c-plane and the energies E$ and 



widths T of the bound and resonant states for the potential V^(r) in several higher partial waves 
(1 < t < 4) . 



t 


fc (fm" 1 ) 


E (MeV) 


T (MeV) 


1 




?2 289054013 


-2.619884138 





1 


1 270932606 


— ?'0 000000043 


0.807634844 


0.000000110 


1 


2.674841953 


- i0.002433228 


3.577386775 


0.013017001 


1 


3 263588553 


— iQ 162882441 


5.312239776 


1.063162540 


1 


3 742982846 


— iQ 564324132 


6.845729429 


4.224511090 


1 


4 225164261 


— iO 999000737 


8.427005280 


8.441884426 


1 


4.696779393 


- il.437816345 


9.996210410 


13.506212364 


1 


5.156711217 


- i\. 873567558 


11.540707589 


19.322893683 


2 




il. 232503483 


-0.759532418 





2 


2 183644493 


— ?'0 000018973 


2.384151637 


0.000082862 


2 


3.052966547 


- ?0.035600651 


4.659668666 


0.217375191 


2 


3.527492840 


- i0.341406094 


6.163323809 


2.408615103 


3 


1.420585762 


- i0.000000016 


1.009031953 


0.000000046 


3 


2.760769155 


- i0.001484354 


3.810922062 


0.008195917 


3 


3.343986231 


- iO. 139396664 


5.581406240 


0.932281051 


3 


3.802701793 


- i0.526397919 


7.091723079 


4.003468620 


4 


2.313665822 


- i0.000006013 


2.676524768 


0.000027824 


4 


3.170315114 


- i0.023354979 


5.025176235 


0.148085287 


4 


3.623738918 


- i0.294896551 


6.522259887 


2.137256217 


4 


4.076486137 


- i0.706384524 


8.059380065 


5.759133441 



TABLE IV. Two lowest Regge trajectories for the potential V^r) . Only those points which 
correspond to the bound and resonant energies are given. 



First Regge trajectory 


Second Regge trajectory 


E (MeV) 


I 


E (MeV) 


I 


-4.571182814 
-2.619884138 
-0.759532418 
1.009031953 
2.676524768 



1 

2 

3.000000000 + i0.000000013 
4.000000000 + i0.000008632 


-0.884280776 
0.807634844 
2.384151637 
3.810922062 
5.025176235 




1.000000000 + i0.000000034 
2.000000000 + i0.000027431 
2.999963480 + i0.003081380 
3.990736977 + iO. 065066318 
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TABLE V. Third and fourth Regge trajectories for the potential V2(r) . These trajectories 
start from the lowest 5-wave resonances. 



Third Regge trajectory 


Fourth Regge trajectory 


E (MeV) 


I 


E (MeV) 


I 


2.252380731 
3.577386775 
4.659668666 
5.581406240 
6.522259887 


-0.000000010 + i0.000041610 
0.999888289 + i0.005370144 
1.976641344 + i0.104936778 
2.805402233 + iO. 392410862 
3.530415810 + i0.760317663 


4.500948186 
5.312239776 
6.163323809 
7.091723079 
8.059380065 


-0.037200122 + i0.130360179 
0.726629120 + i0.451038159 
1.389037872 + iO.838973703 
2.025517176 + il. 253031383 
2.634772674 + il. 671995431 
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FIGURES 




FIG. 1. The zeros of the Jost function for the S'-wave resonances of the potential V\{r) with 
(open circles) and without (filled circles) the Coulomb term. Exact values are given in Table I. 
The dividing line corresponds to the rotation angle 6 = 0.17T. 
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FIG. 2. The zeros of the Jost function corresponding to bound and resonant states generated 
by the potential V^(r) in the S- (open circles) and P (filled circles) partial waves. The filled box 
indicates the D-wave bound state. Exact values are given in Table II and III. The dividing line 
corresponds to the rotation angle = 0.157T. 
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FIG. 3. The five lowest Regge trajectories for the potential V2(r). Exact values of the points 
are given in Table IV and V. Filled circles indicates points which coincide (they are at the same 
place of the £-plane but correspond to different energies) or cannot be distinguished. 
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